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Continuously self-similar solution of spherically symmetric gravitational col- 
lapse of a scalar field is studied to investigate quantum mechanical black hole 
| formation by tunneling in the subcritical case, where, classically, the collapse 

' does not produce a black hole. 
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I. INTRODUCTION 



Recently gravitational collapse of scalar fields have attracted much attention after the 
surprising discovery of critical phenomena by Choptuik The spherically symmetric col- 
lapse of a scalar field has been extensively investigated both analytically and numerically 
The critical solution often displays self-similarity such that there exists a vector field, 
£, whose Lie derivative of the spacetime metric satisfies L^g = 2g. The first known analytic 
solution of this kind was found by Roberts M in the context of counter-examples to cosmic 
censorship, and later discussed by Brady || in connection with critical behavior. Frolov || 
extended Roberts' work in general spacetime dimensions. 

In the work of Roberts-Brady- Frolov |§^| they considered subcritical, critical, and super- 
critical solutions depending on a parameter characterizing the solutions. In the supercritical 
solution a black hole is formed by the collapse, and in the critical case the spacetime is 
asymptotically flat but contains a null, scalar-curvature singularity. In the subcritical case, 
classically, the scalar field collapses, interacts and disperses leaving behind flat spacetime 
without forming a black hole. 

In this paper we reconsider the subcritical case to investigate the quantum mechanical 
black hole formation by tunneling. Although the flux of the scalar field is not intense enough 
to classically form a black hole, it is possible to tunnel through the effective potential barrier 
and to continue collapsing to form a black hole. In the classically forbidden region we 
shall solve the field equations in Euclidean spacetime, whereas Roberts, Brady, and Frolov 
considered only spacetime with Lorentzian signature. We shall show that the black hole 
does form in the subcritical case through the tunneling mechanism. The mass of the black 
hole turns out to be infinity due to the self-similar nature of the solution. But one may view 
the solution as a near-critical solution with a finite mass in the spacetime region of a scale 
much below the relevant large correlation length. 

In the section II, we present the field equations and effective one particle formalism, 
and in the section III we consider solutions in the Lorentzian and Euclidean regions of 
spacetime. There are two kinds of Lorentzian spacetime region : one before tunneling which 
is the one studied by Roberts, Brady, and Frolov. The other is the region after tunneling. 
In the Euclidean region an instanton type solution is obtained. In the section IV we join 
continuously the above three regions to treat black hole formation by quantum tunneling of 
the subcritical collapse of a scalar field. In the last section, we discuss similar problems in 
general dimensions. 

II. EQUATIONS OF SPHERICALLY SYMMETRIC SELF-SIMILAR FIELDS 

We study the minimally coupled scalar field in the four dimensional spacetime whose 
action is 



where K is the trace of second fundamental form of the boundary. Since we are only 
interested in spherically symmetric solutions, we can reduce the action as 




(1) 
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S sph = \J d 2 x^ r 2 [i?( 7 ) + ^{( V r) 2 + 1} - 2( V 0) 2 ], (2) 

where , y a b is the metric in the remaining two-dimensional manifold. It is most easily handled 
in terms of null coordinates such that the metric is expressed as 

ds 2 = -2e 2a dudv + r 2 dn 2 , (3) 

where a and r are functions of both u and v, and dVt 2 is the line element on the unit sphere. 
The Einstein-scalar field equations read 

{r 2 ), uv = -e 2 ° (4.a) 

Off rp p2f7 

2a uv -^yL= 6 —- 2 ^, v (A.b) 

r\ vv - 2cr ;V r, v = -r((p tV ) 2 (4.c) 

r\ uu - 2a )U r >u = -r(<p yU ) 2 (A.d) 

2<P iUV + (kir 2 ) iV (f> jU + (lnr 2 ) )U ^ = 0, (4.e) 

where the last equation is the wave equation for the scalar field 0, and a comma (,) denotes 
a partial derivative. 

In order to find a continuously self-similar solution, we impose the conditions as 

a(u,v) = a(z), r = -up(z), = <f>(z), (5) 
where z is the scale-invariant variable defined as 

The influx of the scalar field is turned on at the advanced time v — 0, so that the spacetime 
is flat to the past of this surface, and the initial conditions are specified by continuity. The 
region of interest is the sector u < 0, v > 0, where we choose signs such that z > 0, p > 0. 
With this choice of the self-similar metric and field, the equations (4.a)-(4.e) become 

p'pz + (pfz - p p = -^e 2CT , (7. a) 



a"z + <j' + P — = -z{i) 2 , (7.6) 
P 



p z 

P -2pa =-p(0) 2 , (7.c) 

p z + 2ap — 2apz = —pz((p ) , (7. a) 

(j)"p + 2(f) p = 0. (7.e) 
Prime denotes the derivative with respect to z. 
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FIGURES 




Figure 1: The height of the potential V(y) is controlled by the value of cq. The regions I, II 
and III are specified. 



With new variables x and y, given by 

x = -hiz, p = y(x)y/z, 
the above equations reduce to, with the metric factor a = 0, 



f = y 2 -2 + ^ 2 
y 

<p = -» , 
y 



(8) 

(9) 
(10) 



where cq is an integration constant, and dot denote derivative with respect to the new 
variable x. 

The equation @ for y formally describes motion of a particle with zero energy in a 
potential 



V(y) 



y 2 



(11) 



We are interested in the production of a black hole by quantum tunneling in the subcrictical 
case, < Co < 1, where the potential, only allows the classical motion of the particle starting 
from y = oo to a minimum value yj. We call this the Lorentzian region I, where Roberts, 
Brady, and Frolov provide explicit solutions in this region. There is another region, we 
call the Lorentzian region II, where classical motion is allowed, for < y < yn- For quantum 
tunneling we need to consider the classically forbidden Euclidean region, y n < y < yj as 
shown in the Fig. 1. 



III. SOLUTION IN THE LORENTZIAN AND THE EUCLIDEAN SECTORS 



(A) SOLUTION IN THE LORENTZIAN REGION I ( Vl < y < oo) 
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In this region where y is restricted as yi < y < oo, the solutions to the equations (|) and 
(0) are 



y = 1 + yl — c o cosh 2x 

1 , { e 2x + T \ 
<f> = — in 



2c V e2x + r ~ 



— lnT 

c 



(12) 
(13) 



where < c < 1, and T = y (1 — c )/(l + c ). We note that the minimum value of y is 
given by 



yj = l + VI -eg, x = 0, (14) 
but there does not occur an apparent horizon as shown by Brady || and Frolov ||. 

r=0 




Minkowski Region 
r=0 



u<0, 



Figure 2: The spacetime diagram is obtained by joining the regions I and II. The presence of 
the region III is indicated by the strip with arrow, whose full structure does not appear in the 
diagram due to its Euclidean nature. 

At the initial boundary (v — 0, x — — oo) we choose the integration constant of <fi such 
that <p(v = 0) = 0, and the stress-energy tensor components are 



T 



uu \ v=0 



T I 

u 5 - L VV\V = 



2 \ 5 - 1 - uv 



T =0 



(15) 



At the final boundary (u = 0,x = +oo) the components of stress-energy tensor are 

1 



T 

uu \u= 



« a (l-cg)' 



T 



n t 

u , -«■ uv 



0. 



(16) 



which implies that the spacetime is flat for u > 0. However we will not consider the full 
region of — oo < x < oo, but limit to the region — oo < x < 0, and connect this to the 
Euclidean tunneling solution with a smooth boundary condition as shown in Fig. 2. 
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(B) SOLUTION IN THE LORENTZIAN REGION II (0 < y < y n ) 

In the region II where y is restricted as < y < yn the solutions to the equations 
and flTt]| ) are 



y 2 = l- VI - eg cosh 2x , (17) 




— lnT, (18) 

Co 



where < T = y (1 — co)/(l + Co) < 1 and < cq < 1. The maximum of y occurs at x = 
such that 



y 2 H = l-Jl-cl x = 0. (19) 



There occurs black hole singularities when y = (r = y^—uv = 0). Here, cosh2x = 
(z + z 1 )/2 = 1/ ^1 — Cq, and so 

z = e 2xo = T or T- 1 . (20) 

The past singularity occurs at z = T, while the future one at z = T -1 . 

Notice that the geometry does not extend to the light cones (u = 0, and v = 0) as the 
past light cone (v — 0, z — 0) is located at the z-value smaller than the allowed limit z = T, 
and the future light cone {u — 0, z — oo) has larger z-value than the allowed limit zq = T~ x . 
There occurs apparent horizons where (v r ) 2 = which implies y + y — 0. At the apparent 
horizon y\ H = Cg/2, therefore 



zh = V 1 "Co or g , (21) 

i — c 

where the first one is the past apparent horizon of the white hole, while the second one is 
the future apparent horizon of the black hole to be formed. 

(C) SOLUTION IN THE EUCLIDEAN REGION (y n < y < yi) 

In this region the geometry is no longer Lorentzian, but becomes Euclidean, which is 
classically forbidden. We consider this region as a tunneling mediated solution, and for this 
purpose we have to change the signature of the geometry and rederive the field equations. 
Fortunately we only need to substitute the scale variable z as a complex number, namely, 

z = e +2- _> z = e ™ ) (22) 

where 9 is the polar angle in the u — v plane. Denoting y = dy/d6, <fi = dcf)/d6 we obtain a 
new set of equations 



2/ 2 = V + 2-4> (23) 

y 

V- (2") 
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where we note that the tunneling potential Vt uri (y) is an upside down flip of the Lorentzian 
potential V(y) in Eq. flU]). 

The solution to the Eqs. fl2|)-(|§ are 

y 2 = 1 + - d cos 29, 0<c <l, (25) 

i J c tan6> 
= — tan j^=^= 

We note that y 2 is periodic in 9 as it should be because it represents a metric in a plane. 
At = 0, y 2 = 1 + \J\ — Cq which is the same value as y 2 in the Lorentzian region I, so 
the Euclidean and the Lorentzian sectors are continuously connected. At 9 = n/2, y 2 = 

1 — y 1 — Cq and this is equal to the y 2 u in the inner region II, here again continuous patching 
of geometry is possible. 



+ 



(26) 



IV. BLACK HOLE FORMATION BY TUNNELING 

In the subcritical case collapsing scalar fields approach only to yi as x approaches to zero 
starting from x = — oo (v — 0), and bounce back to the future light cone (u — 0). However, 
this collapsing scalar field may tunnel through the Euclidean region and appear again with 
the value yn at the matching point with the inner Lorentzian region. For this purpose we 
patch the three solutions to form a continuous complex spacetime as shown in the Fig. 2. 

We recall that in the Lorentzian region I and II, 

ds 2 = -2dudv + r 2 dQ 2 = 2dw 2 - 2w 2 dx 2 + r 2 dQ 2 , (27) 

where 

u = -we' x , v = we x } w > , (28) 

r = wy = y\J —uv . (29) 

The solutions to be patched are 

y 2 = 1 + yl - Cq cosh 2x, - oo < x < 0, (Region I) , (30) 
y 2 = 1 - \Jl-cl cosh2x, < x < oo, (Region II) . (31) 

Notice that at x = these solutions are not continuously matched. What we interpret on 
this mismatch is that, as we can see in the Fig. 2, as x approaches zero the geometry come to 
the classical turning point and it tunnels through the potential barrier along the imaginary 
time while the real time remains at x — 0. In the imaginary time, which is resulted by the 
Wick rotation of the time x, the signature of spacetime becomes Euclidean. Here we use 

ds 2 = 2dw 2 + 2w 2 d9 2 + r W , (32) 

where < 9 < 2n, and x in Eq. (|27j) is substituted by i9. The solution in this region is 
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V 



= 1 + 0" 



Cg cos 20 . 



(33) 



Note that y 2 in the region I and this region continuously joins at x — 0, = 0, their value 
being yj — 1 + J I — Cq. Now the geometry evolves through until approaches 7r/2, where 



y 2 approaches to the value y 2 u = 1 + y 1 — Cq , where it patches to the region II continuously 
at x = 0. 

As shown in Fig. 2, in the spacetime region prior to v — the geometry is flat Minkowski, 
and it evolves along x from x = — oo to x = as the scalar fields collapse, then parts of the 
fields bounce back. But some may tunnel through the potential barrier evolving through 
the Euclidean geometry from 6 = to = tt/2, and the fields reappear at the inner region 
where the spacetime geometry forms an apparent horizon at Xah = — ^ ln(l — Cq), and finally 

collapse to form a black hole at x B h = ^ ln^(l + c )/(l — c )) leaving a singularity (r = 0). 

In order to evaluate the probability to form a black hole by quantum tunneling through 
the Euclidean barrier we follow usual methods of instanton action calculation. We recall 
that the action is 



S sph = \J d'xV^r 2 ^) + ^{(VO 2 + 1} 



2(V# 



udud9((yr) 2 + 1 -r 2 (v0) 2 ) 
where we use the metric, and the coordinates 



^ = ( o -L 2 ) ' u 



-uoe 



ooe 



Since r = uy, we can evaluate (v r ) 2 an d (V0) 2 as 



(yrf = Y'dardtr = ^(y + u^-y) 2 - ^-y 



d 



did' 



;v0) 2 = ^a0) 2 -^ 2 . 



where y — 



The action reads explicitly 



S sp h = ^ / d 1 dx 



-ij 2 + y 2 <p 2 + 2-y 2 + (2q^-y) 2 - V^q^) 2 



(34) 
(35) 

(36) 
(37) 

(38) 



where q = *y and we have dropped an irrelevant total derivative term. 

Now we turn to the particle picture tunneling through the potential barrier as shown in 
the Fig. 1. Defining the momentum densities as 



7T„ 



dC 

dy 



-y, 



dC 



y 2 <i>, 



(39) 



we obtain the effective superspace Hamiltonian density 



H = n y y + 7T, 
2 



- ^ 2y 2 2 



2-l/ 2 + (2^) 2 -y 2 (2^0) 2 



(40) 
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In order to get the equivalent particle picture, one may consistently set ^-\ip) = ^-\ip) = 
by demanding \ip) to be a function only of the zero-mode variables, y = J dqy(q)/ J dq and 
= / dq<p(q)/ J dq. Then diagonalizing tt^ by = \§)e lc ° $ dq<t> with ^|$) = 0, the problem 
is further reduced to 



H\&) 
H = - 



= 

1 d 2 K 
2Kdf + ~2 



y 2 -% 
y 



(41) 



where K = / dq. The effective action of the equivalent particle is then 



S eff = y / d.r 



■J2 

y 



2 + f + ^ 



(42) 



with the Hamiltonian constraint H ~ 0. Notice that the Euclidean version of this action 
gives the equation of motion of the particle given by Eq. (^). The tunneling probability 
can be evaluated within the standard WKB approximation scheme. Inserting the solution 

y 2 = 1 + J I — Cg cos 29 with x = i6, and integrating from # = 0to6' = 7r/2we get 



S, 



eff 




(43) 



where we used K = J d (^J in the last line. It is not possible to use this effective action 
as it diverges unless we introduce a cut-off in the u;— variable. As usual in the self-similar 
solutions this kind of divergences is unavoidable, so we introduce an arbitrary cut-off even 
though this is not allowed in the strict sense of self-similar solution. This divergence is an 
unnatural artifact of self-similarity as the black hole mass also diverges. With such a cut-off 
uj c , we get the probability of tunneling, therefore of black hole production as 



P 



e ~ 2S eff 



3 -^7T(l-Co) 



(44) 



The probability of the bouncing without black hole formation is just 1 — P. We note that 
in the limit of critical case (cq — » 1) we get the expected value P — 1. 

Here the consideration of relative probability to the bouncing case does not remove the 
cut-off dependence contrary to the case of, for example, the free energy in the Schwarzschild- 
Anti-de-Sitter relative to pure Anti-de-Sitter space 0. In order to avoid this divergent action 
naturally we may have to solve the equations with scalar flux of finite duration only, but this 
would require numerical studies. One may interpret the solution as a near-critical solution 
with a finite mass and a finite probability of tunneling in the spacetime region of a scale 
much below the relevant large correlation length. Then the cut-off scale will be naturally 
controlled by the correlation length. 

Brady || defined a local mass function m(u, v) by 
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which agrees with both the ADM and Bondi masses in the appropriate limits. We apply 
this formula to get the black hole mass formed after tunneling (x > 0) as 

cl^uv eg ( -uv \ 2 

1 — Co cosh 2x / 

and, along the apparent horizon (zh = 1/\A — c o? = c o/2); we get 

*U i 

which agrees with the Brady's result ||. In order to compare with the Brady's calculation, 
one should note that his parameters(a, f3) are related to Co as a = (3 = — Cq giving an 
exponent h 

The discrepancy of the critical exponent from the Choptuik's (~ 0.37) may be understood 
as follows. In the near-critical interpretation with a finite but large correlation length scale, 
the apparent horizon will continue to develop in the outside region of the correlation length 
scale, and eventually match with the event horizon if one considers the influx of scalar field 
with a finite duration. The mass viewed in this outside region, then, is expected to change 
appreciably and may conform with the Choptuik's. 



V. DISCUSSION 

Frolov H extended the work of Roberts || and Brady || to the problems in n-dimensions, 
where the action is 

5 = lib/ dTlx ^[R ~ 2(V0) 2 ] , (48) 
which reduces to the spherically symmetric action as 

S sph oc J d 2 x^r"- 2 [i?( 7 ) + (n- 2)(n - 3K 2 (( V r) 2 + 1) - 2( V 0) 2 ] • (49) 
With the self- similarity ansatz we introduce 

d 7 2 = -2e Mz) dudv, r = -up{z), <p = <j>(z) , (50) 
and we get the effective particle motion as 

f = y 2 - 2 + c iy ~^ , (51) 
= c y-^ , (52) 

where 
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Here as in the previous cases we use the definitions, x = |lnz, and p = y/zy(x). The 
solution can be formally presented as 

x = ± [ dy + c 2 , (54) 



which can be performed explicitly in terms of elliptic functions only in the dimension n = 5, 
and n = 6. 

In the Euclidean sector, the solution becomes 

+ c 2 , (55) 



2-y 2 -^ 



y 2(n-3) 

where < 9 < n. The solution y = y{9) is an instanton like solution. We may join the 
above solutions continuously to form a tunneling geometry. Even though the solution need 
not be periodic, it is an interesting observation that the period satisfies the condition, 

<T n <7r (n>5). (56) 



2(n - 2) 



Notice that the period is tt for the four dimensional spacetime, but there are no periodic 
solutions between 5 < n < 9, and yet there exists periodic solutions for n > 10. We do 
not know whether this is somehow connected with supergravity models of dimension ten or 
eleven. 

As a passing remark we mention that the effective potential (|2"5D, V{y) = y 2 + c^/y 2 — 2, 
is a potential of a three dimensional simple harmonic oscillator with an angular momentum 
Cq. It also reminds us the Calogero model || whose Hamiltonian is H = p 2 + X 2 /q 2 . Our 
potential may be viewed as simple harmonic oscillator modified by Calogero potential such 
that solubility might be preserved. 

Finally we should remind readers that our quantum black hole formation process is dif- 
ferent from those black hole pair creations via instanton where gravitational collapse is not 
the driving mechanism 0. Our geometry is more kin to the Hartle-Hawking's no bound- 
ary creation of the Universe flOfl , or Vilenkin's tunneling proposal fll |: our solutions with 



reversed time describes nothing but the creation of the self-similar universe by tunneling, 
although its relevance to cosmology needs to be clarified further. 
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